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ABSTRACT 

This  paper  presents  the  results  of  a  recent  comparison  between  three  microbubble  size  spectrum  measurement 
systems.  These  systems  are  the  light-scattering  bubble  counter,  the  photographic  bubble-imaging  system,  and 
the  acoustic  resonator  array.  Good  agreement  was  formed  among  these  three  systems  over  the  bubble  size  range 
appropriate  for  each  system. 


1.  Introduction 

In  the  past  decade,  the  importance  of  microbubbles 
to  many  physical,  chemical,  and  biological  phenomena 
in  the  near-surface  air-sea  boundary  layers,  both  above 
and  below  the  sea  surface,  has  bwn  well  recognized 
and  documented  (Wu  1987).  It  is  becoming  increas¬ 
ingly  certain  that  bubble  plumes  generated  by  ocean 
wave  breaking  play  a  significant  role  in  underwater 
acoustic  propagation  and  scattering  problems.  The  ef¬ 
fects  of  bubbles  are  seen  in  both  high  ( 10-400  kHz) 
and  low  (50-1000  Hz)  frequency  ranges  through  either 
single-bubble  resonance  or  collected  plume  resonance 
(Prosperetti  1988;  Henyey  1990). 

Unfortunately,  up  to  now  no  single  ideal  technique 
has  been  developed  to  measure  microbubbles  in  either 
the  laboratory  or  in  the  ocean.  An  ideal  system  would 
provide  an  adequate  size  range,  a  laige  sampling  vol¬ 
ume,  a  laige  spatial  area  and  depth  coverage,  and  quick 
response  time,  and  be  amenable  for  automatic  data 
analyses.  The  need  for  such  a  bubble  sensor  system  is 
clear,  and  it  is  still  beyond  our  reach.  At  the  present 
time,  a  combination  of  several  different  types  of  bubble 
sensors  working  in  concert  may  give  an  overall  ac¬ 
ceptable  performance. 

The  Coastal  Systems  Station  (CSS),  since  1984,  and 
Naval  Research  Laboratory  (NRL),  since  1986,  have 
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been  involved  in  the  development  of  sensor  systems 
for  microbubble  measurements  in  the  field.  The  system 
developed  by  CSS  is  the  light-scattering  bubble  counter 
(LSBC),  and  the  two  systems  developed  by  NRL  are 
the  acoustic  resonator  array  (ARA)  and  the  photo¬ 
graphic  bubble-imaging  system  (PBIS).  All  three  sys¬ 
tems  have  been  deployed  independently  to  measure 
wave-generated  bubbles  in  the  sea,  and  they  have  been 
used  in  laboratory  tanks  for  testing  and  calibration  us¬ 
ing  artificially  generated  bubbles. 

It  is  instructive  to  compare  the  performance  of  these 
three  bubble  sensor  systems,  which  operate  on  entirely 
different  optical  and  acoustic  principles,  under  identic^ 
laboratory  conditions.  In  1990,  CSS  and  NRL  jointly 
carried  out  such  a  comparison  test  at  CSS  in  Panama 
City,  Florida.  The  purpose  of  this  paper  is  to  present 
the  experimental  results  of  this  comparison  test. 

In  section  2,  the  operation  principles  of  CSS’s  and 
NRL’s  bubble  sensors  will  be  presented.  In  section  3, 
the  methods  for  the  comparison  test  will  be  described 
with  the  results  presented  in  section  4.  Finally,  a  sum¬ 
mary  and  concluding  remarks  about  these  comparison 
tests  on  bubble  measurement  techniques  will  be  made. 

2.  OpetatioR  principles  of  bubble  sensors 

We  shall  describe  first,  in  detail,  the  operating  prin¬ 
ciples  of  the  three  CSS  r  nd  NRL  bubble  systems,  which 
are  involved  in  the  laboratory  comparison  test.  The 
five  factors  or  criteria  we  shall  use  as  a  base  for  eval¬ 
uating  all  the  bubble  sensors  are 
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FI — bubble  size  range  (radius  in  microns), 

F2 — sampling  volume, 

F3 — spatial  coverage, 

F4 — sampling  time,  and 
F5 — data  analysis  complexity. 

Here,  the  sampling  volume  (F2)  refers  to  the  instan¬ 
taneous  spatial  volume  sampled  by  the  sensor,  while 
the  spatial  coverage  (F3)  refers  to  the  distance  from 
the  sensor  location  that  can  be  measured  by  the  sensor 
remotely.  The  sampling  time  (F4)  is  the  appropriate 
time  period  required  for  obtaining  a  statistically  reliable 
bubble  density  measurement;  thus,  it  is  related  to  ( F2 ) 
and  ( F3 ) .  The  data  analysis  complexity  ( F5 )  is  whether 
or  not  a  bubble  measurements  system  is  amenable  to 
automatic  processing  of  raw  data  to  yield  the  bubble 
density. 

a.  Light-scattering  bubble  counter 

The  light-scattering  bubble  counter  (LSBC)  was 
originally  designed  and  developed  by  C.  S.  Ling,  Cath¬ 
olic  University  of  America,  for  laboratory  investigations 
of  microbubbles  upon  the  inception  of  cavitation  in 
1 982  ( Ling  and  Pao  1 988 ).  It  has  been  further  modified 
into  a  field  system  under  CSS  sponsorship.  In  addition 
to  laboratory  experiments  and  tests,  this  field  system 
has  been  employed  a  few  times  in  sea  operations  (Ling 
and  Pao  1988;  Su  et  al.  1988)  in  conditions  ranging 
from  fairly  calm  seas  in  the  Gulf  of  Mexico  to  violent 
seas  in  the  North  Sea.  Currently,  this  system  is  operated 
by  CSS  personnel. 


This  LSBC  system  is  based  on  dark-held  specular 
reflection  for  detection  of  microbubbles,  and  this  prin¬ 
ciple  is  depicted  in  Fig.  I  [  from  Fig.  3  in  Ling  and  Pao 
( 1988)].  White  light  is  critically  reflected  by  the  gas- 
water  interface  of  a  bubble  that  passes  through  a  sam¬ 
pling  volume  ^  1  of  5  mm  X  6  mm  X  0.4  mm  =  0.0 1 2 
cm^  located  within  the  center  of  an  open-ended  cylinder 
about  3.8  cm  in  diameter.  The  flow  containing  bubbles 
is  sucked  through  this  cylindrical  cavity  by  a  pump. 
The  reflected  light  intensity  is  received  by  a  photo¬ 
multiplier  detector,  at  an  angle  of  125°  with  respect  to 
the  incident  light  and  is  proportional  to  the  square  of 
the  bubble  radius.  The  coefficient  of  proportionality 
can  be  determined  by  the  well-known  relationship  be¬ 
tween  the  bubble  size  and  its  terminal  velocity  in  a 
freely  rising  situation.  This  terminal  velocity  deter¬ 
mination,  during  the  sensor  calibration,  is  facilitated 
by  the  second  sampling  volume  Ai,  which  is  parallel 
toAi  and  separated  by  a  known  distance  (O.S  cm).  The 
simultaneous  recording  of  two  light  intensity  signals 
from  the  two  photomultipliers  are  correlated  to  provide 
the  travel  time  over  the  separation  distance. 

For  better  signal-to-noise  ratio,  the  square  root  of 
the  photomultiplier  voltage  output  is  taken  prior  to 
being  input  to  either  an  FM  recorder  with  a  20-kHz 
bandwidth  or  a  digital  recorder  after  A/D  (analog  to 
digital)  conversion  at  a  rate  of  20  kHz.  Because  of  this 
operation  on  the  data  stream,  the  recorded  signal  is 
directly  proportional  to  the  bubble  radius.  Ling  and 
Pao’s  ( 1988)  Fig.  7  gives  a  calibration  curve  for  the 
LSBC,  showing  that  the  sensor  is  linear  from  a  radius 
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of  about  10  to  about  i  SO  In  other  calibrations, 
the  sensor  is  shown  to  be  useful  for  radius  r  up  to  200 
The  upper  bound  of  r  =  200  is  due  to  the 
thickness  of  the  sampling  volume,  0.4  mm  =  400  fim. 
In  principle,  this  thickness  can  be  increased,  at  the  ex¬ 
pense  of  lowering  the  maximum  bubble  density,  since 
the  sensor  is  allowed  to  have  only  one  bubble  present 
in  the  sampling  volume  at  any  given  time,  otherwise 
causing  an  erroneous  reading.  The  LSBC  can  only  be 
used  in  an  in  situ  mode,  and  its  sampling  time  is  long, 
10  min  or  more,  due  to  its  small  sampling  volume 
(0.012  cm^).  The  technique  is  suited  for  automatic 
data  analyses.  With  a  mean  flow  velocity  within  the 
sampling  volume  of  0.5  m  s  ' ,  the  effective  sampling 
volume  per  10  min  is  a  5  mm  X  6  mm  X  0.5  m  s“‘ 
X  60  s  X  10  min  =  9000  cm^. 

b.  Photographic  bubble  image  system  (PBIS) 

The  photographic  bubble  imaging  system  developed 
by  NRL  consists  of  a  35  mm  camera  in  a  waterproof 
housing  and  three  underwater  strobes  mounted  on  a 
square  stainless  steel  frame,  as  shown  in  Fig.  2.  The 
three  strobes  ( Ikelite )  form  an  equilateral  triangle  with 
each  side  43  cm.  The  strobes  are  located  14  cm  behind 
the  focal  plane  of  the  camera.  This  design  principle 
follows  essentially  that  of  Johnson  and  Cooke  ( 1979). 

The  camera  system  consists  of  the  following  com¬ 
ponents: 

1 )  a  Nikon  F-3  camera  body, 

2)  a  Nikon  55-mm  f/2.8  macro  lens, 

3)  a  Nikon  MD-4  motor  drive 

4)  a  Nikon  MF-4  250  exposure  magazine,  and 

5)  a  Nikon  MF-17  data  back. 

The  lens  is  set  for  a  magnification  of  1:3,  with  an 
f-stop  of  f/ 1 1.  The  focal  plane  is  located  34  cm  from 
the  film  plane  when  the  system  is  submerged,  as  shown 
in  Fig.  3. 


FKj.  2.  PBIS  configuration  used  in  the  comparison  tests. 


To  determine  the  bubble  size  distribution,  each 
photographic  frame  is  viewed  under  a  microscope.  The 
Bausch  and  Lomb  microscope  used  in  the  data  analysis 
has  a  I  OX  objective  and  a  lOX  eyepiece  for  a  total 
magnification  of  1(X)X.  For  bubble  sizing,  a  calibrated 
500-^m  reticle  with  10-;tm  divisions  is  used.  An  X-Y 
scale  ( 1  mm  X  1  mm)  was  fabricated  to  locate  the 
position  of  each  bubble  in  the  frame.  The  bubbles  pro¬ 
duce  a  triangular  pattern  of  dots  that  are  specular  re¬ 
flections  from  the  strobes.  Each  side  of  the  triangle  is 
measured  to  within  5  /im  and  recorded.  These  three 
lengths  are  then  averaged  and  multiplied  by  a  factor 
of  2.25  to  obtain  the  bubble  radius.  The  factor  of  2.25 
is  needed  to  correct  for  the  1:3  image  reduction  and 
other  geometric  conversions.  Figure  4  illustrates  the 
measurement  procedure. 

Since  the  depth  of  field  is  proportional  to  the  bubble 
size,  then 

depth  of  field  (m)  =  0.22  X  0.022  X  radius  (/tm), 

(1) 

with  aperture  at  f/11.  The  bubble  density  for  each 
bubble  size  is  computed  first.  The  total  bubble  density 
is  then  obtained  by  summing  up  the  individual  bubble 
densities.  These  densities  are  then  averaged  over  a  20- 
^m  “bandwidth.” 

Based  on  the  above  relationship,  the  sampling  vol¬ 
ume  per  frame  for  r  =  100  pm  is  about  10  cm  X  15 
cm  X  2.2  cm  =  330  cm^,  which  is  27  500  times  larger 
than  that  of  the  LSBC.  Since  the  strobes  have  a  charge 
time  of  about  5  s  (allowing  12  photographs  per  min¬ 
ute),  the  equivalent  total  sampling  volume  per  10-min 
period  is  330  X  (60  X  10)/5  =  39  600  cm^.  This  is 
about  four  times  larger  than  the  LSBC. 

The  photographic  system  has  to  be  operated  in  situ, 
and  its  sampling  time  is  about  2.5  min,  which  is  equiv¬ 
alent  to  about  25  photographic  frames.  The  bubble 
size  range  for  this  system  is  50-500  ^m  in  radius.  The 
lower  limit  is  caused  by  the  intrinsic  difficulty  in  dis¬ 
cerning  smaller  bubbles  on  a  photographic  negative, 
even  under  a  microscope,  based  both  on  our  own  ex¬ 
perience  and  that  of  Walsh  and  Mulheam  ( 1987).  On 
the  other  hand,  the  upper  limit  is  due  to  the  statistically 
high  uncertainty  associated  with  the  slight  probability 
of  encountering  bubbles  larger  than  this  radius  in  a 
limited  number  of  photographic  frames. 

The  most  serious  drawback  of  the  photographic  sys¬ 
tem,  however,  is  related  to  the  tedious  analyses  of  the 
photographic  negatives.  The  methods  used  by  Johnson 
and  Cooke  ( 1979),  Walsh  and  Mulheam  ( 1987),  and 
our  efforts  presented  here  are  all  manual.  An  automatic 
pattern  recognition  scheme  can  certainly  be  developed 
for  a  digitized  film,  but  this  process  is  still  cumbersome 
due  to  the  great  amount  of  digital  data  per  frame.  At 
present,  we  regard  the  phot(%raphic  system  as  practical 
only  for  calibrating  other  bubble  systems  in  laboratory, 
in  the  bubble  radius  range  of  its  applicability.  For  this 
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special  purpose,  we  can  afford  to  do  a  very  limited 
amount  of  manual  processing,  but  we  definitely  do  not 
recommend  this  system  for  routine  field  operations. 

c.  Acoustic  resonator  array  (ARA ) 

The  acoustic  resonator  array  was  constructed  by 
NRL  and  is  based  on  a  design  originally  developed  by 
Herman  Medwin  and  his  colleague  of  the  Naval  Post¬ 
graduate  School  ( Medwin  and  Breitz  1989;  Breitz  and 
Medwin  1989).  The  system  consists  of  a  capacitive 
transducer  plate  and  a  reflector  plate  held  in  a  parallel 
configuration.  The  transducer  consists  of  a  0.5-mil 
sheet  of  aluminized  mylar  stretched  (not  glued)  over 
an  aluminum  disk  26  cm  in  diameter.  The  aluminum 
disk  is  embedded  in  a  2.5-cm-thick  disk  of  PVC  and 


X 


Fig.  4.  Microbubble  showing  reflection  of  substrobes. 


on  the  front  of  the  mylar  is  a  2-mm  coating  of  Ecco- 
thane,  which  forms  a  waterproof  seal.  The  reflector 
plate  consists  of  a  stainless  steel  disk  6  mm  thick 
mounted  on  PVC.  A  ceramic  disk  hydrophone,  also 
coated  with  Eccothane,  is  embedded  in  the  surface  of 
the  reflector  plate.  The  transducer  and  reflector  are 
separated  by  a  distance  of  24. 1  cm.  Figures  5  and  6 
illustrate  the  construction  of  the  resonator. 

The  operating  principle  of  the  acoustic  resonator  is 
based  on  two  kinds  of  resonances;  the  first  is  the  geo¬ 
metrical  cavity  resonance  induced  by  the  open  space 
between  the  transducer  and  the  reflector,  which  are 
parallel  to  each  other.  When  a  white  noise  is  produced 
by  the  source  transducer  and  reflected  partially  back 
and  forth  between  these  two  plates,  a  series  of  standing 
acoustic  waves  are  created  within  this  open  cavity  be¬ 
tween  the  two  plates.  For  a  separation  of  24. 1  cm  be¬ 
tween  the  two  plates  currently  used  in  our  bubble  sen¬ 
sor,  the  fundamental  frequency  of  the  standing  waves 
is  close  to  3  kHz  and  32  higher  harmonics  also  appear 
in  the  lOO-kHz  bandwidth  of  the  white  noise.  The  sec¬ 
ond  type  of  resonance  is  individual  bubble  resonance 
due  to  excitation  by  acoustic  waves,  which  has  an 
equivalent  cross  section  about  a  thousand  times  latter 
than  the  corresponding  geometrical  cross  section  of  a 
bubble.  For  the  fundamental  frequency  of  3  kHz,  the 
resonant  bubble  radius  is  about  1 200  jtm,  while  for  the 
highest  harmonic  frequency  of  99  kHz,  the  resonant 
bubble  radius  is  about  34  /im. 

If  we  were  to  use  a  400-kHz  bandwidth  white  noise 
instead  of  the  100  kHz  currently  used,  then  the  lower 
limit  of  the  bubble  size  would  be  reduced  to  about  8.5 
ftm.  However,  for  a  frequency  much  higher  than  100 
kHz,  the  acoustic  extinction  (which  includes  both 
acoustic  scattering  and  absorption)  is  no  longer  mainly 
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from  the  individual  bubble  resonance,  but  is  also  from 
the  off-resonance  acoustic  extinction  of  lai^er  bubbles. 
A  much  more  sophisticated  inversion  scheme  such  as 
those  recently  presented  by  Commander  and  Mc¬ 
Donald  ( 1991 )  would  then  be  needed.  On  the  other 
hand,  if  the  separation  between  the  two  plates  is  dou¬ 
bled,  then  the  upper  limit  of  the  bubble  radius  will  be 
doubled  as  well,  to  a  value  of  2.4  mm.  The  price  to 
pay  for  extending  the  range  to  this  larger  bubble  radius 
is  a  poorer  geometric  resonance  between  the  two  plates. 
With  the  above  two  possible  adjustments  and  im¬ 
provements,  the  acoustic  resonator  is  believed  to  be 
capable  of  measuring  bubbles  with  radii  from  about  8 
fim  to  2.5  mm. 

The  complete  block  diagram  of  the  system,  shown 
in  Fig.  7,  was  constructed  and  integrated  by  NRL  for 
use  at  sea.  The  transducer  is  connected  to  a  300-V  dc 
polarizing  battery  and  amplified  white  noise.  The  white 
noise  is  product  from  the  output  section  of  an  HP- 
3561 A  dynamic  analyzer  and  is  amplified  by  an  In¬ 
struments  Incorporate  Model  L2  power  amplifier.  The 
signal  received  by  the  hydrophone  is  fed  into  an  Ithaco 
Model  143E  preamplifier  that  provides  20  dB  of  gain 


and  matches  the  high  impedance  of  the  hydrophone 
to  the  low  impedance  of  the  transmission  cable.  The 
signal  is  then  input  into  the  HP-3561  A  dynamic  ana¬ 
lyzer  where  its  power  spectrum  is  calculated  by  the  fast 
Fourier  transform  { FFT )  and  then  averaged.  The  mea¬ 
surement  time  is  about  equally  divided  for  data  ac¬ 
quisition  and  for  FFT  computation  (Fig.  8).  After  100 
FFTs  are  averaged,  which  takes  about  10  s,  the  dataset 
is  then  downloaded  to  an  HP-9836  or  other  computer 
where  it  is  stored. 

The  conversion  of  the  acoustic  resonator  data  to 
bubble  size  distribution  is  based  on  the  formula  ( Breitz 
and  Medwin  1989) 


Voia)  = 


2irA/ 

culO^Sa  ' 


(2) 


where  jjo  is  the  zeroth-order  bubble  density,  c  is  the 
speed  of  sound  in  water,  a  is  the  extinction  cross  section 
at  resonance,  A /is  the  increase  in  resonance  width,  5 
is  the  damping  constant  at  resonance,  and  a  is  the  bub¬ 
ble  radius. 

The  differences  in  resonance  heights  between  the 
FFT’s  with  bubbles  present  and  with  a  bubble-free  ref¬ 
erence  case  (Fig.  9)  is  then  used  to  compute  the  A/ 
used  in  Eq.  (2).  For  more  details,  see  Breitz  and  Med¬ 
win  (1989).  An  initial  bubble  distribution  no  is  then 
computed  for  each  case.  A  correction  for  the  contri¬ 
bution  of  off-resonance  bubbles  is  then  computed  by 
integrating  the  equation 

agif)  =  -4.3  f  (rv(a)da,  (3) 

VOi 


where  a,  =  34  ftm  and  02  =  1200  fim.  With  n(a)  re¬ 
placed  by  rjo(a)  and  then  substituting  into 


’7(0,)  =  Vo(o,) 


'«g(/,)' 


(4) 


This  produces  the  modified  and  corrected  first-order 
bubble  distribution,  171(a).  This  last  correction  can  be 
carried  to  a  higher  order,  if  desired,  by  repeating  the 
above  procedure. 

The  bubble  radius  range  for  the  current  NRL’s 
acoustic  resonator  array  is  34-1200  /am,  vvdth  a  scat¬ 
tering  volume  about  1250  cm^  as  determined  experi¬ 
mentally  by  Breitz  and  Medwin  ( 1989).  Assuming  that 
the  ARA  near  the  sea  surface  is  under  the  action  of 
orbital  motion  of  ocean  waves  with  2  m  in  wave  height 
and  5  s  in  wave  period,  and  further  assuming  that  the 
mean  flow  rate  passing  through  the  resonator  is  about 
one-quarter  of  the  maximum  orbital  velocity,  the  total 
scattering  volume  over  a  10-min  period  is  about  6 
X  10®  cm^,  which  is  about  600  times  laiger  than  the 
LSBC  and  150  times  larger  than  the  PBIS.  The  ARA 
is  an  in  situ  sensor,  has  a  fast  response  time  of  10  s  or 
less  for  a  reliable  measurement  of  the  bubble  density, 
and  is  amenable  to  fast  automatic  data  processing. 

At  the  present  configuration,  the  NRL’s  ARA  is  op- 
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Fio.  8.  Sample  reference  (bubble-free)  curve  from  the  acoustic  resonator  0-100  kHz. 


erated  in  the  field  by  being  tethered  to  a  research  vessel 
by  a  cable  about  1 50  m  long.  But  the  system  can  be 
easily  converted  into  an  autonomous  system  with  bat¬ 
tery  power  for  about  100  h  of  continuous  operation. 

3.  Comparison  tests 

The  comparison  tests  among  the  CSS’s  LSBC  and 
NRL’s  PBIS  and  ARA  were  conducted  at  CSS,  Panama 
City,  using  two  freshwater  test  tanks  in  1990. 


The  first  test  was  between  ARA  and  PBIS  and  the 
second  test  was  between  ARA  and  LSBC.  In  both  tests, 
we  employed  three  methods  for  generating  bubble 
plumes  (see  Figs.  lOa-c).  The  first  method  (Fig.  10a) 
used  electrolysis  with  the  bubble-producing  wire 
wrapped  around  a  30-cm-long  rod  that  was  made  to 
sweep  over  80°  to  generate  wider  and  more  homoge¬ 
neous  bubble  plumes.  The  second  method  produced 
bubbles  by  using  a  series  of  15  hypodermic  needles 


FK3.  9.  Sample  reference  and  bubble  curve  from  the  acoustic  resonator  0-200  kHz. 
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Fig.  10.  Three  methods  employed  for  bubble  generation;  (a)  electrolysis  wire  wrapped  on  a  sweeping  rod, 
(b)  compressed  air  passed  through  needles  on  a  sweeping  tube,  and  (c)  underwater  jet  air  entrainment. 


mounted  on  a  copper  tube,  which  is  connected  to  a 
compressed  air  tank  (Fig.  10b).  This  tube  assembly 
also  sweeps  in  the  same  way  as  the  first  method.  The 
third  method  used  an  underwater  jet  for  entraining  air 
into  the  water  in  such  a  way  (Fig.  10c)  that  a  large 
bubble  plume  is  created.  The  bubble  size  range  created 
by  the  underwater  jet  is  wide,  from  tens  of  microns  to 
several  millimeters.  The  bubble  density  may  be  con¬ 
trolled  to  a  certain  degree  by  varying  the  flow  rate  of 
the  input  tap  water  and  by  adjusting  the  jet  head  lo¬ 
cation  and  direction  relative  to  the  sampling  volume 
of  bubble  sensors. 

After  analyzing  and  comparing  of  the  bubble  den¬ 
sities  measured  by  these  three  bubble  sensor  systems 
using  these  three  different  methods  for  providing  bub¬ 
ble  plumes,  we  found  that  the  third  method  of  under¬ 
water  jet  is  the  most  satisfactory  for  providing  more 
continuous  bubble  density  over  the  range  from  30  to 
1200  fim  than  the  first  two  methods.  For  this  reason. 


all  the  discussions  on  bubble  tests  shall  be  restricted  to 
those  being  employed  by  the  underwater  jet  only. 

During  the  first  comparison  test  between  ARA  and 
PBIS,  the  two  systems  were  hung  from  two  supporting 
beams  on  the  top  of  the  test  tank,  with  PBIS  placed 
above  the  ARA,  as  shown  in  Fig.  1 1.  A  bubble  plume 
was  generated  about  65  cm  below  the  ARA.  Thus,  the 
same  bubble  plume  rose  through  the  sampling  volumes 
of  the  two  bubble  sensors  that  were  centered  about  3.0 
and  3.4  m,  respectively,  below  the  water  surface.  Data 
were  collected  simultaneously  by  both  sensors  under 
this  condition  during  the  test. 

During  the  second  comparison  test  between  LSBC 
and  ARA,  again  both  sensors  were  hung  in  the  water 
tank,  as  shown  in  Fig.  12,  in  such  a  way  that  the  flow 
intake  opening  of  the  LSBC  was  locat^  about  5  cm 
above  the  top  edge  of  the  ARA.  Again,  the  same  bubble 
plumes  pas^  through  the  sampling  volumes  of  the 
two  bubble  detectors  that  were  centered  about  1 .6  and 
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Fig.  1 1.  Configuration  of  NRL's  ARA  and  PBIS 
used  in  comparison  test. 


1 .95  m,  respectively,  below  the  water  surface.  The  water 
in  the  entire  test  tank  was  made  to  be  fully  saturated 
with  air  prior  to  the  test  so  that  the  effects  of  natural 
degasing  of  air  bubbles  into  the  surrounding  water  was 
minimized  during  their  free  rise  from  the  first  sampling 
volume  at  the  bottom  to  the  second  sampling  volume 
at  the  top.  The  pressure  differential  effect  on  the  bubble 
size  is  very  small  (about  1.5%  of  the  bubble  radius) 
because  the  separation  between  the  two  scattering  vol¬ 
umes  is  only  about  35  cm. 

Since,  in  both  series  of  comparison  tests,  the  centers 
of  the  two  scattering  volumes  of  the  pair  of  bubble 
sensors  were  at  different  water  depths,  any  single  bubble 
will  change  its  radius  traveling  from  one  volume  to  the 
next.  The  ambient  pressure  effects  on  bubble  size  have 
been  corrected. 

Each  of  the  two  comparison  tests  were  repeated  ten 
times  by  the  first  and  second  methods,  and  five  times 
with  varying  bubble  densities  by  the  third  method  of 
bubble  generation. 

The  test  data  collected  were  processed  according  to 
the  principles  and  formula  as  described  in  section  2  and 
finally  expressed  in  the  bubble  density  in  number  of  bub¬ 
bles  per  cubic  meter  per  micron  range  of  the  bubble  radius 
versus  bubble  radius  in  the  customary  log-log  plot  for 
accommodating  the  large  dynamic  range  of  the  density 
over  a  radius  range  of  10-1000  ^m. 

4.  Results  of  the  comparison  tests 

We  shall  present  first  results  of  the  test  between  the 
acoustic  resonator  array  (ARA)  and  the  photographic 
bubble  imaging  system  (PBIS).  Figures  13,  14,  and  15 
are  three  (out  of  five)  cases  of  this  comparison.  The 
other  two  cases  show  similar  results.  In  each  plot  we 
show  the  bubble  densities  measured  by  the  two  systems 
under  the  same  conditions.  Overall  good  agreement 


between  these  two  bubble  densities  over  the  radius 
range  of  80-400  fun  is  observed.  The  larger  fluctuation 
seen  in  the  bubble  density  by  the  PBIS  is  believed  to 
be  due  to  an  insufficient  number  of  photographic 
frames  taken  and  analyzed,  that  is.  statistical  fluctua¬ 
tion  due  to  undersampling.  Furthermore,  the  bubble 
densities  at  radii  less  than  80  //m  is  somewhat  lower 
for  PBIS.  This  is  attributed  to  the  reduced  sampling 
volume,  caused  by  the  smaller  depth  of  field,  and  to 
the  relative  insensitivity  of  the  film  for  smaller  bubbles 
(Walsh  and  Mulhearn  1987).  On  the  other  end  of  the 
bubble  radius  range,  the  bubble  density  by  the  PBIS  is 
also  lower  than  that  by  the  ARA.  This  might  be  also 
caused  by  statistical  undersampling,  since  there  are  only 
a  few  large  bubbles  in  the  frames. 

Figure  16  shows  three  successive  bubble  densities 
each  averaged  over  10  min  by  the  ARA  using  the  un¬ 
derwater  jet  for  bubble  generations.  Very  good  consis¬ 
tency  among  them  is  observed  over  the  radius  range 
from  34  to  600  ^m.  This  clearly  shows  the  steadiness 
of  the  bubble  plume  generation. 

Next,  we  present  the  results  of  comparison  tests  be¬ 
tween  the  ARA  and  the  light-scattering  bubble  counter 


Fig.  12.  Configuration  of  ARA  and  LSBC 
used  in  comparison  test. 
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Fig.  13.  Comparison  of  bubble  densities 
measured  by  PBIS  and  ARA. 


( LSBC ) .  Three  ( out  of  five )  bubble  densities  obtained 
by  the  two  different  sensors  using  three  variations  in 
flow  rates  of  the  underwater  jet  are  shown  in  Figs.  17, 
18,  and  19.  The  other  two  cases  show  similar  results. 
One  can  see  clearly  from  these  figures  that  the  bubble 
densities  obtained  by  these  two  methods  agree  remark- 


Fta.  14.  Comparison  of  bubble  densities 
measured  by  PBIS  and  ARA. 


RadiMiMa] 

FKj  .  I S.  Comparison  of  bubble  densities 
measured  by  PBIS  and  ARA. 

ably  well  in  the  radius  range  of  34-200  (im,  which  is 
currently  the  applicable  radius  range  of  the  LSBC,  while 
the  applicable  range  of  ARA  extends  to  1200  ^m. 

Among  the  three  bubble  sensors  compared — that  is, 
ARA,  LSBC  and  PBIS — we  have  reasons  to  think  that 
the  LSBC  is  the  most  accurate  sensor  in  the  bubble 
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Fig.  17.  Comparison  of  bubble  densities  generated  by  water  jet, 
measured  by  LSBC  and  ARA. 


radii  range  from  30  to  200  jtm,  given  the  condition 
that  bubble  plumes  are  statistically  steady  and  that  the 
measurement  time  is  sufficiently  long.  These  reasons 
are  as  follows:  First,  the  bubble  density  is  obtained  in 


Fig.  18.  Comparison  of  bubble  densities  generated  by  water  jet 
measured  by  LSBC  and  ARA. 


Fig.  19.  Comparison  of  bubble  densities  generated  by  water  jet. 
measured  by  LSBC  and  ARA. 


the  LSBC  by  a  simple  counting  of  individual  bubbles 
in  a  definite  and  well-defined  sample  volume.  Second, 
the  determination  of  the  bubble  radius  for  those  bub¬ 
bles  are  based  on  the  quite  accurate  terminal  velocity 
relationship  with  respect  to  bubble  size,  confirmed  in 
the  laboratory,  resulting  in  a  very  linear  calibration 
curve.  Third,  the  ARA  uses  some  physical  approxi¬ 
mations  in  deriving  the  inversion  formula  and  has  a 
less  well-defined  sample  volume  that  is  not  uniformly 
active  with  respect  to  acoustic  extinction  and  that  must 
contain  many  bubbles.  Taking  these  factors  into  con¬ 
sideration,  we  conclude  that  the  bubble  densities  mea¬ 
sured  by  the  LSBC  in  its  applicable  range  ( 30-200  ^m) 
are  the  most  reliable  reference  data  that  we  have.  The 
very  close  agreement  between  the  bubble  densities 
shown  in  Figs.  1 7- 1 9  by  the  ARA  and  LSBC,  in  effect, 
proves  the  nearly  equal  accuracy  of  the  ARA  method, 
far  beyond  the  more  limited  supporting  evidence  orig¬ 
inally  given  by  Breitz  and  Medwin  ( 1989).  The  smooth 
and  regular  extension  of  the  bubble  density  beyond  the 
r  =  200  nm  in  these  figures  further  strongly  suggests 
that  the  accuracy  of  the  ARA  may  be  taken  up  to  r 

=  1200  j«m. 

5.  Summary  and  conclusions 

We  (NRL  and  CSS)  have  jointly  conducted  two  se¬ 
ries  of  laboratory-controlled  tests  in  order  to  compare 
the  performances  of  a  photographic  bubble  imaging 
system  (developed  at  NRL),  a  light-scattering  bubble 
counter  (developed  at  CSS),  and  an  acoustic  resonator 
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array  (developed  at  NRL).  Three  kinds  of  bubble  gen¬ 
eration  methods — electrolysis,  hypodermic  needle  ar¬ 
ray,  and  underwater  jet  air  entrainment — were  em¬ 
ployed  to  provide  bubble  plumes  in  the  test  tank  com¬ 
mon  to  any  pair  of  two  bubble  sensors  under  test.  We 
found,  however,  that  the  underwater  jet  air  entrainment 
provided  the  most  suitable  bubble  plumes  for  our 
comparison  tests.  As  such,  we  presented  only  the  results 
of  using  the  bubble  generation  method  in  this  paper. 
Bubble  densities  were  obtained  from  all  three  bubble 
sensors  under  the  same  bubble  plumes.  Comparisons 
of  these  bubble  densities  shows  that  over  the  bubble 
radius  range  from  about  75  to  400  ftm,  the  results  from 
photographic  bubble  imaging  system  agree  well  with 
the  acoustic  resonator  array,  but  that  beyond  the  above 
range,  the  former  deteriorated  rapidly  due  to  photo¬ 
graphic  film  insensitivity  and  statistical  undersampling. 
This  comparison  further  shows  that,  over  the  radius 
range  of  34-2(X)  ^m,  the  light-scattering  bubble  counter 
agrees  remarkably  well  with  the  acoustic  resonator  ar¬ 
ray.  Since  among  these  three  bubble  sensors,  the 
acoustic  resonator  array  has  the  largest  sampling  vol¬ 
ume,  fastest  response  time,  and  widest  radius  range  (at 
present  from  34  to  1 200  /nm  and  can  be  extended  to  8 
/im  to  2.5  mm)  and  is  the  easiest  to  operate  and  perform 
automatic  data  analyses,  it  is  obviously  the  most  suit¬ 
able  method  for  bubble  density  measurements  in  the 
oceanic  upper  layer.  As  a  matter  of  fact,  we  have  de¬ 
ployed  three  times,  in  1991-92,  the  NRL’s  ARA  in 
the  oceans  for  bubble  density  measurements  under 
moderate  to  storm  conditions  (wind  speed  10-17 
ms"').  Results  of  these  field  measurements  are  being 
analyzed  and  are  to  be  published  in  the  near  future. 
The  bubble  measurement  closest  to  the  surface  is  about 
25  cm,  the  same  as  reported  by  Medwin  and  Breitz 
( 1989).  However,  in  some  laboratory  situations  with 
more  confined  test  space  and  steady  bubble  plumes, 
the  light-scattering  bubble  counter  is  unquestionably 
more  suitable. 


It  is  too  often  that  bubble  sensors  are  developed  and 
used  without  having  any  comparison  with  existing 
sensors  under  carefully  controlled  situations.  In  con¬ 
ducting  these  comparison  tests,  both  agencies  (NRL 
and  CSS)  learned  valuable  lessons  concerning  the  three 
sensors  tested. 
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